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Photodisintegration of the deuteron above tt threshold is studied in a coupled channel approach 
including A'^A and vrd channels with consideration of pion retardation in potentials and exchange 
currents. A much improved description of total and differential cross sections in the energy region 
between tt threshold and 400-450 MeV is achieved. With respect to polarization observables, the 
description of the linear photon asymmetry and the proton polarization remains problematic. 



I. INTRODUCTION 



Recently we have constructed in Q, henceforth denoted as I, a realistic hadronic interaction model in the two- 
nucleon system which takes into account meson retardation completely. This model is based on meson, nucleon 
and A degrees of freedom, and it is able to describe NN scattering quite satisfactorily also above pion threshold. 
^ However, for technical reasons it only allows configurations where maximal one meson is present explicitly at a time. 
■ Within this restriction, full meson retardation in the potentials has been incorporated in a consistent manner. The 
, motivation for this work was twofold: First of all, from the general principle of relativity any interaction has to 
00 ■ be retarded taking into account the finite velocity of its propagation. With respect to the NN interaction, this is 
particularly evident for energies above pion threshold, where an exchanged pion can become an onshell particle, so 
• that the retarded meson propagator becomes singular describing the coupling of the NN channel to the now open 
\ ttNN channel. Secondly, there is strong evidence in deuteron photodisintegration above tt threshold, in particular 
in the A resonance region, that neglected retardation is one major cause for the failure to describe experimental 
data by quite sophisticated theoretical treatments which, however, are based on static interactions and corresponding 
static exchange currents. Indeed, first results have demonstrated clearly the considerable improvement of the theory 
if retardation is considered |§,D. 
(-H ^ In the present work, we have studied deuteron photodisintegration in greater detail with respect to the influences of 
the various ingredients of the theoretical framework developed in I. In Sect. |l| we discuss the general structure of the 
e.m. current operator, which consists of the baryonic components, for example, the one-body spin, convection, and 
isobar currents, as well as of mechanisms where mesons, especially pions, are created or absorbed. The latter ones are 
^ ' a natural consequence of the incorporation of mesons as explicit degrees of freedom. As is well known, the internal 
absorption of a photon by a nucleon creating an intermediate A isobar is of particular importance in the A resonance 



(3/2) 

the M^_^i multipole in pion photoproduction on the nucleon. In this way, all parameters of the th eory are fixed in 
advance before studying electromagnetic reactions in the two-nucleon system. In the following Sect. |^ the resulting 
effective current, acting in the baryonic space only, is derived. Special emphasis is devoted to the structure of the 
retarded MEC which has, in contrast to its static counterpart, a singular structure above tt threshold describing the 
coupling of the NN to the ttNN channel. In Sect. ^ the question of gauge invariance is discussed. The results for 
total and differential cross sections as well as for selected polarization observables are presented in Sect. Finally, 



region. Therefore, Sect. [II is devoted to a careful determination of the corresponding 77VA coupling by studying 



Sect. VII contains a summary and an outlook. 



II. THE ELECTROMAGNETIC CURRENT OPERATOR 



The transition amplitude Tfi, which describes the absorption of a photon by a hadronic system making a transition 
from an initial state \pi,i) to a final state \pf,f), is given in Coulomb gauge by the Fourier component of the e.m. 
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current matrix element Jfi{k) 

S^^\pf - P. - k)Tf, = -e{k, A) • Jf,{k) , (1) 
where the polarization vector of the incoming photon with momentum k and helicity A is denoted by e(fc, A). 

A. The general structure of the current 

We begin with a brief discussion of the formal structure of the e.m. current. As has been described in detail in I, 
the model Hilbert space of the two-nucleon system consists of three orthogonal subspaces 

^[2] = ® n^l^ ® H^x^ (2) 

[2] [2] [2] 

containing either two bare nucleons (H^), one nucleon and one A (7i^ ), or two nucleons and one meson {Ti-x )- 
Accordingly, the current operator can be decomposed into various diagonal and non diagonal components. With the 
help of the projection operators P^, Pa, and Px onto the above subspaces, respectively (see I), the current J^(fc) 
can be written as a symbolic 3x3 matrix 



( 4A(fc) 4x(^) \ 

V "^XN^^^ -^XA^^) "^XxW J 



(3) 



Now, we will discuss the various components in detail. We distinguish between pure "baryonic" currents, "one-meson 
production/annihilation" currents, and the remaining currents, and divide the first one further into one- and two-body 
operators, labeled by a superscript "[1]" and "[2]", respectively. The diagonal baryonic currents J'^^{k) and J^^{k) 
are represented diagrammatically in Fig. 0. The two-body parts comprise effective exchange currents, like for example, 
heavy meson exchange currents, which are not generated explicitly via the elementary meson-baryon vertices and the 
meson production/annhilation currents. The transition current J^^(fc), also represented in Fig. describes the e.m. 
excitation of a A resonance including two-body contributions. 

The other transition current consists of three components, namely jx]^\ ^xff^^ -^xt? ^' "^here the super- 
script "(0)" or "(1)" indicates the order with respect to the meson-nucleon coupling constant. Of these, the one-meson 

'^'^ and i^l")'^ 



production currents j^]^ ^-^^^ ^xN ^ shown in Fig. |^. Note that the current Jxa ^^'^^ present because it is set 



equal to zero (see the discussion in Sect. |[VD[) . In detail, j^j^^ denotes the contact current which is related to the 



KroU-Rudermann term of pion photoproduction, and jx^'^ vertex-current, which is generated by the hadronic 

vertex form factors. We include in J'^j^{k) in general only the pion as most important contribution to the retarded 
exchange currents with one exception, namely those parts which are necessary for generating the retarded ^np/uj 
MEG. This current is of some relevance in deuteron photodisintegration and therefore is taken into account. The 
remaining current jx]^^ which is a pure mesonic transition current, is shown in diagram (a) of Fig. ^together with 
the two one-body parts of the diagonal component J^jj-(fc), describing the coupling of a photon to a nucleon (meson) 
with a spectating meson-nucleon (nucleon- nucleon) state (see diagrams (b) and (c) , respectively) . These contributions 

are denoted as jxxi^) ^^"^ Jxx(^)' respectively. The current jx]^ gives an important contribution to the retarded 
MEC which, however, violates the one-meson-approximation, being another pathology of this approximation which 
we had alluded to in L But neglecting it, would lead to a severe violation of gauge invariance for the retarded tt MEC. 
Therefore, it has to be included perturbatively in the construction of the corresponding effective exchange operator. 

B. ExpHcit form of the current components 

Most of the currents are obtained by the canonical method of minimal coupling, and the corresponding matrix 
elements of J are derived in a straightforward manner and need not be given here explicitly. The nucleon one-body 
current, appearing in jj^^^ and jxxW^ needs some special consideration. Since we start from bare nucleons, the 

corresponding current jj^'j^ is given by a Dirac current with vanishing anomalous magnetic moment k. However, as 
discussed in detail in I, the bare nucleon becomes dressed by meson-nucleon loop contributions (see Appendix A of I). 
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Consequently, the dressed nucleon acquires a nonvanishing anomalous magnetic moment by the presence of this meson 
cloud (see Fig. However, it is not in agreement with the experimental value as needed for a realistic treatment of 
photonuclear reactions. Therefore, we substitute the dressed nucleon current of the model by the onshell current with 
the experimental value for k. The corresponding nonrelativistic current is denoted by i"g'^['(fc) and contains the usual 
convection and spin part. In addition, we also take into account the relativistic spin-orbit current as most important 
relativistic contribution of leading order. This relativistic contribution is quite important in photodisintegration of 
the deuteron, even at low energies 

It is obvious that the above substitution of the physical onshell nucleon current creates formally several inconsis- 
tencies. For example, the neglect of the explicit evaluation of the loop diagrams leads to a violation of unitarity if the 
nucleon current operator is evaluated in the two- nucleon system where the initial and/or final nucleon can become 
offshell. In this case, some of the loops implicitly present could become singular above pion threshold and, therefore, 
would lead to an imaginary part describing real meson production. 

Now we will consider the remaining current contributions. The matrix elements of the vertex current j^^x ^ 
given by [D 

(7V(p')|pi^l(fc)k(<f,/i)^(p)) = 0, (4) 
{Nip')\4'x^ {k)H<l, t^)Nm = Siq+ p + k - p') ^ [e, r,] {2q + k) ^a ■ {q + k) ^-(g') ' fc) ') ^ (5) 

where e = e (1 + T3)/2, and a pion state with momentum q and isospin projection /i is denoted by |7r((f, /i)}. As next, 

we identify the nucleon current contribution in the presence of a spectator meson, Jxxik), with the normal nucleon 
current 

{niq,^^)NifM^ik)\niq\f,')N{p)) = V (2^)' Mq) 5{<i' ~ ?) {N {p'Wlai {p)) ■ (6) 

The current components for the "/Trp/uj currents cannot be obtained by minimal substitution, but from the usual 
interaction Lagrangian 

= ^pe^.A. (a--^^ - a^vi--) (a^'/'^ - aV") 0. with xe{p,u;}. (?) 

The electromagnetic coupling constants g-yirx, as determined from the decay of the vector meson into a pion and a 
photon, arc |^,^ 

g^^p = 0.41 , g^^^ = 2.02 . (8) 

Now, we will turn to the currents in which a A isobar is involved. For the direct electromagnetic exitation 
7 -|- TV — > A we take the usual nonrelativistic form which contains a magnetic dipole (Ml) and an electric quadrupole 
{E2) transition. The analysis of pion photoproduction on the nucleon indicates that the E2 transition is largely 
suppressed compared to the dominant Ml transition. Therefore, we neglect the E2 contribution and thus take as 



current |10 



qOAN 

(A(pA)|jAJv(fc)l^(pjv)) = e jv. '^(PA - Pjv - fe) JJ} i^AN x kjN , (9) 

ZlVIx 



with 



kjN^k ^——JLp^^ where M^^^ = 1232MeV. (10) 

The corresponding jNA charge, which is associated with the current contribution proportional to pA, is of minor 
importance and therefore neglected in this work. The Ml strength is determined by fitting the M-^^ ' multipole 
amplitude in pion photoproduction on the nucleon which will be the topic of the next section. 

At the end of this section we will fix the remaining current contributions J^^(fc'), J^J^(fc) and j]^j^(fc). Due 
to the absence of A7V components in the deuteron wave function, the current J^^{k) cannot contribute to the 
electromagnetic deuteron break-up. jj^J^(fc) is identified with the p MEC J^'^^^{k), which is the most important two- 
body current besides the tt MEC. Note, that first of all in jj^|^(fc) as well as in jj^j^(fc) bare couplings occur, which 
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have to be renormalized (see Eq. ^ in Sect. IVC). In our approach, j'^^^^ik) is taken in the strictly nonrelativistic 
form [0. One should note that, due to the neglect of retardation in J ^P\k), this MEG is not fully consistent 
with the p exchange of the Elster-Potential which serves as basic NN interaction in this work (see I). Similarly, the 



second contribution jj^j^(fc) is identified with the static pionic A MEG J^'{k) in the nonrelativistic limit In 
photodisintegration, this current contribution is almost negligible. 



III. THE Mf^'^'' MULTIPOLE OF PION PHOTOPRODUCTION 

In order to fix the coupling constant G^J^^ for the Ml excitation of the A resonance, we now will consider the 

(3/2) 

multipole of pion photoproduction on the nucleon. The production amplitude is given by the sum of a Born 
and a resonance part (see Fig. 

i.^(z)=tf°™(z)+t?-(z). (11) 
The resonance part is given in compact form by 

t^-(z)=«X„AgA(^)JAivW, (12) 

where denotes the dressed A propagator as introduced in I. The effective, energy dependent 7iVA current Jan{z) 
contains a direct and a rescattering part, i.e. 

N[i] 

where go denotes the nucleon propagator, and N^i^ a renormalization constant in the one-nucleon sector as defined in 
the Appendix A of I. 

The dominant contributions to the Born amplitude to the m}^^^' amplitude are displayed in Fig. |^. They consist 
of the time ordered pion pole diagrams and the crossed nucleon pole diagram. Following the procedure in [|l3|, less 
important mechanisms like u exchange will be taken into account effectively by a parameter b which is close to unity. 
Thus we use explicitly as Born amplitude 

(^(q»7V(-<f)|if;<^7„(z)|7(fc)^(-fc)) - -6^[e,rt]e.(2g- k)ta.{q~ k) F^{{q - k)^) 



2M' 
1 



2uo{q- k) 



1 



z ~ k~ uj{q~ k ) — e'pj'iq) z — oj{q) — i^{q — k) — e^{k) 



, y-n eTle-{2q + k)ia-q+ {e + k)Tle-a xka-q 
4A'^Ar z - k - uj{q) - e'^{q + k) 



where the form factor Ft, and the pion-nucleon coupling constant are given by the underlying NN interaction 
(here the Elster-potential, see Table I in I). The parameter h and the renormalized 77VA coupling constant = 
G^'^^/7V[i] are fitted to the experimental data of Arndt et al. jl^ (solution SM97K) yielding 

G^P^ = -3.91 and b 0.84. (15) 

We would like to emphasize that Watson's theorem jl^ is exactly fulfilled below two-pion threshold. In view of small 
but significant differences of the experimental Pas-phases determined from TriV scattering versus pion photoproduction, 
we have determined G^H and h as well as the hadronic parameters f^T^j^, ^AttN, and by a simultaneous fit of 

the m[^^^ multipole and the P33 partial wave of ttN scattering (see Figs. and ^). 

As next, we would like to introduce the effective strength G^^ (z, k) which is defined by the matrix element of the 
effective current J^pfiz) in the A rest frame 

^ ^ - ^ G^^iz k) 

(A(pA - 0)\j^^{z,k)\N{p^)) =eT^f,^oKPN + k) ' la^n x k , (16) 

where for initial onshell nucleons z and k are related by z ~ W + ie with W — e^{k) + k. In the later discussion 
of the results this coupling will be referred to as (efFl). Because of the occurrence of intermediate pion-nucleon 
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loops in the rescattering amplitude (see Fig. ||), G^,^ (z, k) becomes complex above tt threshold. Modulus han{W, k) 
and phase ^{W,k) of the effective coupling 

Gif = W + ze,k) = J1an{W, k) e^*(^''=' (17) 

are shown by the full curves in Fig. ^. Obviously, the energy dependence of the modulus is rather weak. 

Finally, we would like to compare briefly the effective ^N/S. coupling of Wilhelm et al. 1 12 1^ . Similar to 0,|l8| , 
the onshell matrix element of the resonant amplitude together with the Born terms are fitted directly to experimental 
data without explicitly evaluating the rescattering amplitude. The parameters h and -Pn- are set equal to one, and the 
energy dependent ansatz is used 

MAiv(M^,fc) =Mo + M2 f— ) and$(T^,fc) = , (18) 

where W = e^{q) + w^((f ) and q denotes the momentum of the outgoing pion. The free parameters are fitted to the 
data of Berends and Donnachie [|l9|. This coupling, denoted henceforth by Gf/^ (eff2), is represented by the dotted 
curves in Fig. ^. Note that the modulus of this coupling is considerably smaller at low energies than the one of our 



approach, where the rescattering amplitude is explicitly evaluated. As will be seen in Sect. VI, this feature turns out 
to be one of the reasons for the considerable underestimation of the total cross section of deuteron photodisintegration 
in IT^ra. 



IV. THE EFFECTIVE CURRENT OPERATOR IN DEUTERON PHOTODISINTEGRATION 

After having fixed all current components of (^, we will now construct an effective current operator, which acts in 
pure hadronic space, by eliminating all explicit meson d.o.f. It is defined by 

{f\pj^^^{z,k)p\t) = {f\j'^m). (19) 

The initial and final states are given by the effective deuteron and outgoing NN scattering states, respectively. Their 
explicit forms are given in (55) and (47) through (49) of I. We will evaluate the effective current in the antilab system, 
where the final NN scattering state is at rest and the deuteron has total momentum —k. In principle, one would then 
need a boost contribution for the moving deuteron which, however, can safely be neglected according to the findings 
in H]. 

Inserting the various components of the current into the rhs of (|l9[), it turns out that the effective current operator 
can be split into a nucleonic and a resonant part with superscripts "N" and "A" , respectively, which in turn can be 
divided into one- and two-body terms denoted by superscripts "[1]" and "[2]", respectively (see Fig. |l^) 

j^ffiz, k) = jf^w ^(z, k) + ^(z, k) + j^}f^(z, k) + ^(z, %) . (20) 

We would like to remark, that in this evaluation we have included perturbatively the contributions from the one- 
meson-approximation violating diagram (a) of Fig. |[ In this context, we have to introduce into the corresponding 
expressions below projection operators PxiX2 on ^ state consisting of two nucleons and two mesons of type xi and X2- 



A. The nucleonic one-body contribution J^j-^^ ^ {z,k) 

According to the discussion in subsection |ll B| , this component is represented by the physical onshell current in- 
cluding the relativistic spin-orbit part 

^-\NN; p,a)\PJ^ff^{z,k)P\d) = —^-HnN- p,a)\P^J^ ^ iCanhk)+r'{^,k))P^\d) , (21) 

Nd R{z) ' 

where a characterizes additional quantum numbers, and the normalization constant TV^ is given in Eq. (56) of I. It 
differs from unity in retarded approaches only (e.g., N^ — 0.992 for the Elster-potential). 
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B. The resonant one-body contribution J^'^'^(z, fc) 
The component J^j'^(z, k) is given essentially by the effective jNA current 
^~\NN;p,a)\PAjtff''{^,k)P^\d) = ^^^-\NN;p,a)\PA ^ z - eN{Pj},k)) R^z ~ k)Pj^\d) . (22) 

Note that similar to the Elster-Potential, we consider only pion-nucleon loops in R^{z — k). The effective current 
J^j^(z) has the same functional form as given in except that G^,f[^ is replaced be the effective excitation strength 
^Mi i^sub + ie, fc^jv), where 

->2 

Wsub (W,pa) = W- Mjv - (23) 

denotes the invariant mass of the t:N subsystem which we have evaluated for the spectator onshcU choice. In the 
actual evaluations in the two-nucleon system, G^/^ is calculated in the rest frame of the A, similar to the treatment 

(3/2) 

in the multipole. The difference with respect to the one-nucleon sector is that k^N and Wsub are not related 

by an onshell relation, i.e., e^lkj^) + fc-yjv 7^ Wsub- Because of the explicit evaluation of the pion-nucleon loop in 
the rescattering Born amplitude (cut A in diagram (c) of Fig. ^), this feature does not lead to any complications in 
practice. On the other hand, in the treatment of Wilhelm et al. such an offshell extrapolation of GfH is not 
possible. Therefore, they have used an onshell prescription 

r^Nrycr , ,v A- .A ^ / '^mHEa) for Ea = Wsub > + Mn , 

Gmi {Wsub + k^N) ^ I Gfj^MN + m^) otherwise . ^^^> 

It turns out that this onshell treatment for is a very good approximation, at least for photodisintegration, 

compared to the exact calculation, where the arguments k^N and Wsub in (H^) are treated independently. 



C. The nucleonic meson exchange current J^'J' ''(z, fc) 
The matrix element of J^|^'^(z,fc) has the following structure 

^-HNN-p^a)\P^J^}f^{z,k)P^\d) = ^i~\NN-p,a)\Pf,J^}f^{z,k)^^^^P^\d), (25) 

where J^j^^^{z, k) consists of a static p MEG, a retarded tt MEC, and a retarded 77rp/ti;-current jC**'*)^^^, fc) 

J^/f''iz,k) = J('')''(fc) + J("^'^(z,fc) + J(*""^^(z,fc). (26) 

A new feature of retardation is the fact, that above pion threshold the retarded currents can be split into a part 
containing a pole and a regular part (nonpole) 

ji.yi,.ss^,^^j^^jMms^^ (27) 

The pole part describes the reabsorption of a real pion by a nucleon. Furthermore, the retarded tt MEC splits into a 
contact (c), pion- in-flight (/), vertex (v) and recoil contribution (r) 

jW A?) = j('^)/"^(z, fc ) + j('^)/-'^^(z. A?) + j('^)/"''(z, fc ) + jW/'^^(z, fc ) . 

(28) 

A graphical illustration of these currents is given in Fig. |ll|. In addition, because of the nonvanishing ird interaction 
Vxxj currents like those depicted in Fig. |2| would have to be taken into account. However, a detailed analysis has 
shown that these contributions can be neglected in deuteron photodisintegration. This approximation is realized by 
the substitution 
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G^(z) = (z - H^^xx - V^x) ' Go(z) - (z - Ho.xx)'' (29) 

in the above expressions (see Eq. (49) of I). 

Explicitly, one obtains for the various components 

4oie^'i^^k)^ v%^^{j)Goiz)j^^l'^{k,i), (30) 

ji:Jj;oiei^, k)= E (^^ go(^ - fc)4,^(j-) , (31) 



jM/ft^ 



{z,k) = E (4x,(j)Go(^)jt^l(fc)Go(z-fc)4^^(*) 

+4^„(j)Go(z)^;lx,J0Go(z)j(,%(fc)) , (32) 

Jnonpilei^, k) = E ^fxj^) G,{z - fc^^^ (z)Go (z - fc)4,^(j) > (33) 

^S/''"(->fc)= E 4;f.(j-)Go(^)jS'(fc,^), (34) 

4:itre(-. k)- E - k)v%^^{,) , (35) 



j(7r)/rp 



(^.fc)= E (4x.(jOGo(^)j;Taz^(fc,o4.^(0 + 4xJOGo(^)Car(^.^)4,^(j)) , (36) 

4:itre(-.fc)=o, (37) 

(z,fc) = E (4x,(0Go(z)jtJl(fc)Go(z-fc)4^^(j) 

+4xJj)Go(^)4,x,.(*)Go(^)ji"^'^^(fc)) ^ , (38) 



j{diss) ^ 
'^pole 



ij = l,2:j#i 



+4x.(*)Go(^) A"il(fc) Go(z - fc)4^^(j) 
+4^^ (z)Go(z)4^^^^ (j)Go(z) J^°i''^(fc) 

(fc)Go(z-fc)4^^^^(j)Go(z-fc)4^^(i)) (39) 

In the present evaluation, we have used a nonrelativistic reduction of the irN vertex for the sake of simplicity. The 
bare couplings have to be renormalized, of course, with the renornialization operator Z"^^ in (25). Therefore, it is 
useful to introduce renormalized MECs according to 

J{p)/M/id^ss)^(^^^^ ^ _J_jip)/M/{diss)p^^ J:^J_ (40) 
^[2] ^[2] 

SO that the relevant matrix element fP5l) has the structure 



(-)(AriV;p,a|PJ,^fl^(z,fc)P|d) = (iViV; p, ajP^^ ( j'") ^(^) 

+k'(z) (j(")^(z, fc) + j(*^^)^(z, fc )) i?"(z - fc))p^M') , (41) 

where the dressing operator K" is neglected in the p MEC due to its minor importance. 

Summarizing the various contributions to the nucleonic meson exchange currents, we include as retarded currents 
the TT and ^irp/oj MEC, whereas of the heavier mesons only the p MEC will be retained in the static limit. All 
the other remaining cr, 5, u, and r] MECs are completely neglected in view of their negligible contributions found 
previously in ||2l|. Moreover, due to the use of Siegert operators |2|] and the much larger mass of the heavy mesons, 
the role of the corresponding MEC is largely suppressed. 
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D. The resonant meson exchange current J^'J' ''(z, fc) 

In analogy to the nucleonic MEC, we construct the retarded A MEC which is depicted in Figs. |l^ and The 
currents represented in Fig. |l^ are the contact, meson-in-flight, vertex, and dissociation contributions to the A MEC. 
In the static limit, they are already taken into account via the static A MEC (apart from the dissociation current) 
so that the problem of double counting arises. In view of the fact that this static A MEC leads to a very small 
contribution in deuteron photodisintegration, it is resonable to assume that retardation will not change this fact 
qualitatively. Therefore, we take the diagrams of Fig. |l3| in the static limit as already contained in the static A MEC 
and consequently set J^'^(A;) equal to zero. 

The recoil contribution of Fig. ^ on the other hand, is not contained in jj^j^(fc). Therefore, it is explicitly taken 
into account. Introducing again rcnormalized MECs, we obtain for the matrix element of the A MEC 

^-\NN;p,a\PAJ^ff''{z,k)Pj^\d) = -^J-^ {NN; p, alP^^J^^]^^ (k) + ji^''^'' ' {z,k)R^z ~ k)^^ (42) 
where the quantities J^^f^ik) and j'^^^^^{z, k) are defined in analogy to (|4C|). 



E. The effective current operator in the static limit 



Until now, we have only considered the retarded form of the effective current. In the static case, considerable 
simplifications occur. Whereas the one-body contributions remain unaltered (apart from the dressing factors which 
are, of course, not present in static calculations because in this case only physical instead of bare nucleons appear), 
we have to substitute the retarded meson- A^A^ propagators in Eqs. ( pO| ) through ( p9[ ) by their static limits, i.e. 

Go(z) and Go(z-fc) — > -7^ , (43) 



where hx denotes the kinetic energy of a meson x. Consequently, the pole structure of the MEC is lost. The meson- 
nucleon vertices in the corresponding expressions are those used in the static Bonn-OBEPR potential ||2^ so that the 
sum of the pion-in- flight, vertex and contact MEC fulfils current conservation with the vr exchange part V'qbbp/?, 
the Bonn-OBEPR potential 



{j'stat (k) 



+ Jstaiik) + J Stat (k) ) — 1^0S_EPi? : Preal{k) 



(44) 



Note however, that a static recoil contribution does not appear, because in the static limit the wave function renor- 
malization current cancels exactly the recoil current |2^,^^. Likewise, the recoil contribution in the A MEC (Fig. [l^ ) 
has to be neglected, too. 



F. Inconsistencies in static approaches 

At the end of the discussion of the effective current operator, we would like to point out some inconsistencies in 
static approaches, namely the fact that the A current and the tt MEC are not independent from each other. From 
Fig. ^ it becomes evident that the Born terms of the M^J^"^^ amplitude are related to a part of the tt MEC in the 
two-nucleon system, namely to a part of the pion-in- flight and the recoil current. A static treatment of the MEC leads 
therefore to serious inconsistencies. For example, the crossed nucleon pole graph is not contained in the static tt MEC 
because the respective recoil MEC is not present at all in static approaches due to its cancellation against the wave 



function renormalization current. How critical this inconsistency is has been demonstrated by Wilhclm et al. ^2 16 

'Ml 



in determinin g an other effective 7iVA coupling, denoted henceforth by (effS), in a fit to the data where the Born 



amplitude in (|ll|) is set equal to zero so that it is contained effectively in this coupling. At least below the resonance 
position this "modified" coupling is considerably larger as the original coupling (compare the dashed and dotted 
curves in Fig. In view of this rather large difference between the two couplings, Wilhelm et al. had already 
suspected that retardation in the MEC may become important in the A region. Using the modified parameters, a 



much improved description of the total cross section was achieved (see Sect. VI), though the differential cross section 
was still poorly described in the A region. Moreover, it is obvious that the modified coupling is conceptually not very 
satisfying because the neglect of the Born amplitude leads formally to a vanishing rescattering amplitude, so that 
the resulting effective 7iVA coupling is identical to the bare one and therefore energy independent, in contrast to the 
present parametrization in Eq. ([iT 
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V. THE QUESTION OF GAUGE INVARIANCE 



Now, we would like to discuss the question whether the effective current J^j^{z,k) is gauge invariant, i.e., is 
consistent with the hadronic interaction. As first, we will formulate the continuity equation for the effective operators. 
From the continuity equation for the full current operator 



{f\([H,f\k)\ - k-J{k)) |z)=0 
one obtains straightforwardly the relation 

{f\Pk- Iffizf, k) P\i) - {f\P (H,ff{zf)p,ff{zf,k) - peffizf, k)Heff{z,)) P\t) , 



(45) 



(46) 



where Zi/ f = Ei/ f ±ie, and the interaction part Veff{z) of the effective Hamiltonian consists of a disconnected and a 
connected part 



with 



v,ffiz)^vf;i{z) + v,Yf{z), 



V:ff{z) = V^^;^ + [V^M^)V^ 



XP\ 



Note that we have neglected Vxx i'^ K™/ view of the discussion in Sect. IV C . 

Splitting the hadronic interaction and the current into one- and two-body parts, one obtains 

{f\Pk-j}}\izf,k)P \z) = {f\P [{HoPP + Kffizf)) pW^(z/, k) - pW/z/, fc) (i/opp + P N) 



(47) 

(48) 
(49) 

(50) 
(51) 



These equations serve as starting point for the discussion of gauge invariance. We will restrict ourselves to the pure 
nucleonic currents and interactions acting solely in for the sake of simplicity. Concering the one-body part in 
(pO|), it is clear from the very beginning that in retarded calculations, where V^"^*! is not zero, gauge invariance cannot 

exactly be fulfilled. Note that we do not evaluate explicitly the loops in Fig. ^ which enter into J^^^^^iz, k) as has 
been discussed in subsection [IB. On the other hand, the corresponding loop in Vf^j has to be evaluated due to the 
requirement of unitarity of the hadronic interaction. Therefore, Eq. (^0|) is only valid in the onshell case, i.e., for 



(52) 



where Pi/j is given hy E^/f = 2e^{pi/f). Note that in this case Vfjj is exactly zero. But in static approaches, 
Eq. ( pO| ) can always be fulfilled because of Vfjj: = 0. 

We now turn to the discussion of the two-body part in (|5l|). According to (p6|), the effective current operator 
jjv[2] fi^^^ consists of the tt and p MEG and the jTrp/uj currents. With respect to the p MFC, we refer to the discussion 
at the end of subsection 



II B 



Concerning the ^np/uj contribution, one should notice that the corresponding charge 
operator vanishes due to the use of nonrelativistic vertex structures, whereas the remaining spatial part j('''^*^(z, k) 
is purely transverse 



k- J^'^'''\z,k) = 0. 



(53) 



Consequently, (|5l|) is fulfilled with respect to the dissociation current. The remaining retarded tt MFC satisfies 

(Z, ^Zf ~k) 



k ■ j^-\zf,k) = v:ff-{zf)pruk) - preaiik)Krriz,) + kj^-^°{zf,k) + A{zf, k)+d (^^) 

with the retarded tt exchange potential 



(54) 
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L J con 

- note that in ( [55| ) y^jo"'"^' occurs due to the use of the nonrelativistic pion-nucleon vertex in the MEC -, and the 
auxihary quantity 



- {Ef - H,^^) (pf^Jk) Go{z k)v\-:jf{i)G,{z - k)vl^^{j)) } . (56) 



In general, only between plane waves the matrixelement of A{z, k) vanishes, which fact indicates that in a retarded 
approach additional MECs of at least fourth order in the pion-nucleon coupling are necessary for ensuring gauge 
invariance. Examples for such contributions are depicted in Fig. ^ They are quite complicated and beyond the 
scope of the present work. In the static limit such MECs are not present because they either cancel out each other or 
cancel against corresponding wave function renormalization currents so that (|5^) is satisfied. 

In summary, gauge invariance is fulfilled up to second order in the meson- nucleon coupling constant, whereas 
violations appear in higher order. However, we hope that by the use of the Siegert decomposition of the electric 
multipoles p3] at least a part of these neglected higher order contributions are taken into account implicitly. 



VI. RESULTS FOR DEUTERON PHOTODISINTEGRATION 



We now will discuss the results for deuteron photodisintegration. The formalism for calculating the various observ- 
ables may be found in the review [ p2[ . As currents wc consider all currents discussed before, i.e., the one-body spin and 
convection current as well as the spin-orbit current and tt and p MEC. The parameters of the coupling constants are 
chosen consistently with the underlying NN interction. We have included all multipoles up to L < 4. For the electric 
transitions we have used the Siegert operators which incorporate a large fraction of MEC contributions. They 
include besides the usual one-body density also two-body charge terms as described in Sect. IV C . The remaining 
MEC contributions beyond the Siegert operators are evaluated explicitly. In view of the fact, that retardation is most 
important at short distances, we have considered retardation in the contact, vertex, meson- in- flight and dissociation 
MECs only for multipoles up to L < 2, while for L > 2 the static expressions are used, which is sufficient as we have 
checked in detail. 



A. Retardation effects in a pure nucleonic model 

We will begin the discussion by studying first the influence of retardation in a pure nucleonic model, where the A 
d.o.f. as well as the nd channel are switched off so that only a pure NN interaction and the corresponding currents 
are taken into account. For the sake of simplicity, the dissociation current is also switched off. In detail, we consider 
four different types of calculations which differ with respect to static or retarded potentials and MECs, comprising 
contact, vertex and meson-in-flight currents, and with respect to the recoil MEC. We will denote them by N(potential, 
current, recoil current), and list the four cases which we consider in Table |. As static potential we use the Bonn 
OBEPR model, while for the retarded interaction we take the Elster model. 

It is obvious that due to the neglect of explicit A d.o.f. no realistic description is possible above about kiat = 100 
MeV. But here we are interested mainly in the relative sensitivity of the cross sections with respect to the various static 
and retarded approaches. For this reason we show in Fig. ^ only the ratios of the total cross sections with respect 
to N(stat,stat,0). One readily notices that retardation in the hadronic interaction, i.e., going from N(stat,stat,0) to 
N(ret,stat,0), leads to a drastic reduction of the cross section which amounts to about 25 percent at kiab — 240 MeV. 
On the other hand, retardation in the meson-in-flight, contact and vertex MECs is almost negligible. It turns out 
that retardation is important only in the lowest multipole transitions to the ^5*0 and ^Pq partial waves, which are 
of minor importance in the total cross section, except directly at threshold. In general, the influence of retardation 
in the MEC decreases rapidly with increasing multipole order. On the other hand, the recoil current contribution 
turns out to be very important as is evident comparing N(ret,ret,0) with N(ret,ret,l) in Fig. |l^. It is of particular 
importance in the contributions of the ^Z?2 and partial waves. 

In summary, it turns out that below tt threshold, the static and the full retarded calculations predict the same 
cross section within a few percent. This result is at variance with the previous work of Schmitt et al. |^^, in which 
retardation in potential and current had been considered using a p/Mif expansion below tt threshold. A sizeable 
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decrease of the cross section was found in retarded calculations compared to a static approach. In contrast to the 
present work, also the wave function renormalization current was included but no renormalization of the continuum 
wave functions was applied. Only the bound deuteron wave function had been renormalized globally by a constant. 
In a correct approach, however, the wave function renormalization current should have been left out as well as the 
renormalization of the deuteron state. This erroneous treatment lead to the mentioned decrease of the cross section. 

In the energy region above tt threshold, the equivalence between static and retarded frameworks breaks down, the 
difference increasing with increasing energy as can be seen in Fig. |l^. This is not surprising, because above tt threshold 
the pion can become an onshell particle which fact is only taken into account by the retarded ttNN propagators, 
whereas for a static framework no coupling between the NN and the open t:NN channel is possible, i.e., retarded 
and static approaches exhibit a different behaviour. 



B. Retardation effects in the A resonance region 

Now we will turn to the energy region of the A resonance, applying our full model which includes explicitly A 
degrees of freedom. In detail, we will consider four different static and five different retarded approaches whose 
nomenclature is explained in Tables || and III . The corresponding hadronic interaction models are described in detail 



in I and are therefore not repeated in this paper. We have restricted the full inclusion of A d.o.f. to partial waves 
with total angular momentum j < 3, and use for > 4 - similar to fl^ - the impulse approximation [27| in which 
the A'^A components are given in first order perturbation theory by 

PaIA^A^; p,a)^-^ = G^,^\z) V^^°Ij^{z)Pj^\NN; p,a)^-^ , (57) 

where the nucleonic component \NN; p,a)^ ^ is obtained in pure nucleonic space from a realistic NN potential. 
As first we consider the results of the static MEG cases in Figs. |l^ and 19. The model CC(statl), which has already 



been used by Wilhelm et al. |16|, shows the already mentioned significant underestimation of the total cross section 
by about 30 % at kiab — 260 MeV. Moreover, above kiab — 300 MeV the differential cross sections develop a dip 
structure around 90° in clear contradiction to the experimental data. Due to the slightly stronger ^N/S. coupling in 
the model CC(stat2) - compare the couplings G^H (eff 1) and G^.j^ (eff 2) in Fig. ^~ this model leads to a noticeable 
enhancement of <Ttot which is, however, not sufficient to bridge the gap to experiment. The additional improvements, 
in the approach CC(stat3), especially the incorporation of the 77rp/w currents, lead to a further enhancement of the 
cross section for kiab ~ 340 MeV. Compared to the recent Mainz data , the theoretical underestimation of atot 

amounts to not more than 10 % at kiab = 260 MeV but increases to about 40 % at kiab = 440 MeV. Moreover, 
the dip problem in the differential cross section has not been resolved. The last static model CC(stat4) shows, as 
has already been discussed in ||l^, a rather good description of the total cross section in the peak region due to the 
substantially stronger 7A^A coupling, but again the dip problem as well as the underestimation of the cross section 



at higher energies remains. Moreover, as has been stressed in Sect. IV F, this approach is quite unsatisfactory from a 
conceptual point of view. 



The results for the retarded current models are shown in Figs. 20 and |21|. The effect of introducing retardation 
in the hadronic interaction but keeping the currents still static can be seen by comparing CC(stat3) with CC(ret4). 
One readily notes a drastic decrease of the cross section. This fact has been noted already for the pure nucleonic 
case discussed in the previous subsection. Switching on in addition retardation in the currents (transition CC(ret4) 
— » CC(ret2)) leads, however, to a drastic increase which overcompensates the reduction of the cross section by the 
retardation in the hadronic interaction. It turns out that within about 10 % an overall satisfactory description of 
the total cross section over the whole energy region is achieved. Concerning the differential cross sections, apart 
from energies below the resonance maximum, the description of the experimental data is almost quantitative, too. 
Especially the above mentioned dip structure occurring in the static treatment vanishes completely. This is mainly 
due to the additional consideration of the A^- A- mass difference in ^^"^ ^aA' ^^-'^ "^^^ to the recoil current 

contributions. On the other hand, similar to the pure nucleonic case, retardation in the pion-in-fiight, vertex, contact 
as well as in the dissociation currents is negligible in the total and differential cross sections. With respect to the 
infiuence of the ird channel, it turns out that this mechanism leads to a noticeable enhancement below about kiab — 
300 MeV. This result is at variance with the results of Tanabe and Ohta |^9| who found a decrease of atot when they 
incorporated the 7rd channel. On the other hand, the inffuence of the ird channel on the channel is qualitatively 
the same than in our approach. Thus, the reason for this difference is unclear. 

As next topic, we discuss in Fig. ^ the influence of the parameter aAWp describing the additional pA^A coupling in 
Eq. (76) of I. Comparing the results of CC(retl), CC(ret2) and CC(ret3), which differ by the value of aAWp, but which 
are equivalent concerning the description of the channel in A'^A'^ scattering, it turns out that the value of a^Np is 
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of considerable importance for energies above about kiab = 380 MeV. But similar to NN scattering, an optimal value 
for ttAATp could not be determined because of the rather large differences between the different experimental data sets 
at higher energies ||2^,^. 

Finally, we discuss a few polarization observables. At first, we will consider the linear photon asymmetry S for 
which new experimental data is available |3^j3^. In Fig. |2^, the results for the approaches CC(stat3), CC(ret4) 
and CC(ret2) are depicted. The static approach CC(stat3) yields quite a good description for energies up to about 
kiab = 300 MeV. For higher energies, however, an oscillatory structure appears which vanishes completely if retardation 
in the hadronic interaction is switched on (CC{stat3) CC(ret4)). Moreover, E decreases drastically with increasing 
photon energy. The additional inclusion of retardation in the currents (CC(ret4) CC(ret2)) leads to an overall 
reduction of E except for the highest energy for which a significant change of the shape appears. Compared to the 
experimental data one finds a satisfactory overall description of the qualitative features. But one notes also some 
significant deviations around 90°, where theory predicts an asymmetry considerably smaller than experiment. This 
is particularly significant at kiab — 220 MeV in view of the small error bars, while for higher energies the larger errors 
make the deviations less significant. It is interesting to note that for kiab — 220 MeV also the differential cross section 
exhibits around 90° a sizeable deviation of the theory by about 10 %. 

Finally, we show in Figs. ^ and ^ the proton and neutron polarization, respectively. At kiab — 240 MeV, the 
different approaches yield very similar results which are also in good agreement with the experimental data. At 
kiab = 400 MeV hadronic retardation results in a strong overall decrease of the polarization for both neutron and 
proton (transition CC(stat3) CC(ret4)). Switching on the current retardation leads however to a strong polarization 
increase in the forward and backward direction while around 90° the polarization is much less affected. This feature is 
clearly at variance with the few proton polarization data. Due to the fact that in general polarization observables are 
more sensitive to small reaction mechanisms, one reason for this failure may be the fact that apart from the dominant 
channel several other partial waves (for example the P- and the ^F^ waves) are described only fairly well within 
the present model (see I for further details). Thus, one has to await future improvements of the hadronic interaction 
model in order to see whether this failure could be resolved. Certainly, more data of higher accuracy are needed in 
addition. 



VII. SUMMARY AND OUTLOOK 



In this paper we have presented the electromagnetic part of a model which is able to incorporate complete retardation 
in the two-body n MEC and thus is suited for the study of electromagnetic reactions in the two-nucleon sector above 
pion threshold up to about 500 MeV excitation energy. The electromagnetic current of the hadronic system comprises 
the nucleonic one-body currents, the A excitation and meson exchange currents. The electromagnetic excitation of the 
A resonance has been fixed by fitting the M-^_^^ multipolc of pion photoproduction. Within the present framework, 
it is possible to construct the necessary retarded tt MEC as an effective two-body current operator. In this context 
we would like to emphasize, that we do not use any wave function renormalization [p^ , because this procedure is 
applicable in retarded models only for energies below the tt threshold. Consequently, recoil contributions besides 
the meson-in-flight, vertex and contact MEC have been taken into account in addition. The retarded tt MEC is, 
however, only partially consistent to the corresponding retarded NN interaction. This is mainly due to the neglect of 
- technically very complicated - MEC of at least fourth order in the ttNN coupling constant which are not present 
in static approaches. On the other side, due to the use of the Siegert operators for the electric multipoles, such 
explicitly neglected mechanisms are at least partially taken into account. We furthermore would like to emphasize 
that in the present model all free parameters are fixed in advance before studying electromagnetic reactions on the 
two-nucleon system. In order to test the quality of the developed framework, we have studied photodisintegration 
of the deuteron, where until now serious problems existed in the theoretical description above tt threshold. While 
below TT threshold, say up to about 100 MeV, both frameworks, static and retarded, are equivalent, this equivalence, 
however, breaks down when approaching tt threshold and even more so above, e.g., the total cross section in the 
retarded framework is considerably enhanced compared to the one in the static framework. Within about 10 %, 
we obtain a satisfacory description of the experimental data over the whole energy region of the A resonance. This 
improvement is the combined result of various additional independent mechanisms beyond the ones already considered 
in JTot . These additional ingredients comprise besides the very important retardation effects the dissociation currents, 
the TTd channel and the explicit evaluation of the rescattering amplitude in fixing the 77VA coupling. Furthermore, 
also the theoretical description of the differential cross section is considerably improved. Thus, for the first time a 
model is available which is able to describe the total as well as the differential cross section up to about kiab — 450 
MeV. However, some problems in desribing the linear photon asymmetry as well as the proton polarization at higher 
energies are still present. 
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TABLE I. Nomenclature for the pure nucleonic approaches of Sect. 



VIA. 



notation 


potential 


MEC 


recoil current 


N(stat,stat,0) 


OBEPR 


static 


not included 


N(ret,stat,0) 


Elster 


static 


not included 


N(ret,ret,0) 


Elster 


retarded 


not included 


N (ret, ret, 1) 


Elster 


retarded 


included 



TABLE II 

hadronic model follows Ref. 



Nomenclature for the static models of the coupled channel approach of Sect. VI B 
Q. The recoil current is always not included. 



The nomenclature for the 



notation 


hadronic model 


MEC -yNA coupling j'*""' " 


CC(statl) 
CC(stat2) 
CC(stat3) 
CC(stat4) 


CC(statl,7r) 
CC(stat2,7r) 
CC(stat,7r,p,0) 
CC(statl,7r) 


static Gf/l{eS2) not included 
static G^/l (eSl) not included 
static G^^(effl) included 
static G^/^(eff3) not included 


TABLE III. Nomenclature for the various retardeii models of the coupled channel approach of Sect. VI B. The nomenclature 
for the hadronic model follows Ref. In all cases, the 'yNA coupling G^/f (eff 1) has been used. Except for CC(ret4), the 
retarded recoil current is always included. 


notation 


hadronic model 


MEC remarks 


CC(retl) 
CC(ret2) 
CC(ret3) 
CC(ret4) 
CC(ret5) 


CC(ret,7r,/5,-l) 
CC(ret,7r,p,0) 
CC(ret,7r,p,l) 
CC(ret,7r,p,0) 
CC(ret,7r,p,0) 


retarded 
retarded 
retarded 

static Jf'^'/'', J^^''^'' not included 
retarded nd channel not included 
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FIG. 1. Diagrammatic representation of the baryonic currents. A shaded ellipse symbolizes a two-body exchange current. 





(a) 



(b) 



FIG. 2. Diagrammatic representation of the meson production currents J^^ 



(a) contact current j^S'J^ and (b) vertex current 



(a) 



(b) 



(c) 



FIG. 3. Diagrammatic representation of the two- meson production current j^^]^ (a) and the current components J^^ : (b) 
nucleon current jxxW' i^) meson current jxxW- 
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(d) 
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(f) 



+ 
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(g) 



(h) 



(i) 



FIG. 4. Diagrammatic representation of the separate contributions to the physical one-nucleon current, represented in con- 
trast to the bare nucleon current by a filled circle: (a) bare rmcleon current, (b) and (c) contact current, (d) - (f) meson-in-flight 
current, (g) bare nucleon current with meson in flight, (h) and (i) vertex current. 
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FIG. 5. Diagrammatic representation of the resonant amplitude of pion photoproduction (a) consisting of a direct amplitude 
(b), and a rescattering contribution (c). 




FIG. 6. Diagrammatic representation of the Born amplitude t^"'^{z) in the Mi+(3/2)-multipole amplitude. The separate 
diagrams are denoted as follows: time ordered pion pole terms (a) and (b) and crossed nuclcon pole term (c). 
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1100 1150 1200 1250 1300 1350 
W [MeV] 

FIG. 7. Experimental Pas-scattering phase ^33 of ttN scattering (triangles from [Q, solution SM95) and photoproduction 

phase 5 (3/2) (open circles from [|l4|, solution SM97K). The solid curve represents our fit using the values of Eq. (67) in I, in 
1+ 

which both data sets are considered with equal weight. 
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FIG. 8. Fit of the real and the imaginary part of the m[^^-' multipole amplitude to the experimental data from solutii 
SM97K. The dotted curve represents the contribution of the Born amplitude. 
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FIG. 9. Energy dependence of the effective yNA coupling for the onshell-case W — e^ik) + k. The left panel shows the 
modulus MajvC^i^) ^-nd the right panel the phase ^{W,k). The full curve represents the coupling G^/^(effl) in (|l3|), and the 
dotted curve the one of Wilhelm G^j f (eff2 ) (see @). The dashed curve shows the "modified" coupling of Wilhelm et al. 
Gmi (6ff3) (see the discussion in Sect. IV F). 
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FIG. 10. Diagrammatic representation of the separate contributions to the effective current operator J^^^(2, k) (see Eq. (^)): 



(a) Jf™''(^,fc), (b) Jf'!l^(^,fc), (c) J^ff^iz,k), and (d) jf'^l^(^,fc). 
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FIG. 11. Diagrammatic representation of the retarded n (contact (c), pion-in-flight (f), vertex (v), recoil (r)) and "fnp/cj 
MEG (diss). The arrows indicate the static limit. 
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FIG. 13. Diagrammatic representation of the retarded A MEG: (a) meson-in-flight, (b) contact, (c) vertex, and (d) dissoci- 
ation contributions (a; € {p, a;}). 
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FIG. 14. Diagrammatic representation of the retarded A recoil MEG J^^^^^{z,k). 
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FIG. 15. Correspondence of Born terms of pion photoproduction to tt MEG. 




FIG. 16. Diagrammatic representation of retarded MEGs which are of fourth order in the pion-nucleon coupling constant. 
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FIG. 17. Ratios of total cross sections of d{'y,N)N using pure nucleonic models with respect to crtot(N(stat,stat,0)): dotted 
curve: N(ret,stat,0), dashed curve: N(ret,ret,0), full curve: N(ret,ret,l). 
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FIG. 18. Total cross section of deuteron photodisintegration in the various static approaches listed in Table O: full curve 

CC(statl), dashed curve CC(stat2), dash-dotted curve CC(stat3), and dotted curve CC(stat4). Experimental data from |2^] 

(.), (□) and g (o). 
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FIG. 19. Differential cross sections of deuteron photodisintegration in the static approaches for various photon energies kiab 
as function of the cm. proton angle 9p. Notation of the curves and experimental data as in Fig. |l^. 
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FIG. 20. Total cross section of deuteron photodisintegration in various static and retarded approaches. Notation of the 
curves (see Tables |l^ and ^for the nomenclature): dash-dotted CC(stat3), dashed curve CC(ret4), full curve CC(ret2), dotted 
curve CC(ret5). Experimental data as in Fig. [isl 
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FIG. 22. Differential cross sections of deuteron photodisintegration for various values of qajvp and various photon energies. 
Notation of tlie curves: dashed curve CC(retl), full curve CC(ret2), dotted curve CC(ret3). Notation of the experimental data 
as in Fig. [18. 
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FIG. 23. Linear photon asymmetry E of deuteron photodisintegration. Notation of the curves: dash-dotted CC(stat3), 
dashed CC(ret4), full CC(ret2). Experimental data from jil) (o), Q (•) and jas) (□). 
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FIG. 24. Proton Polarization Py{p) of deuteron photodisintegration. Notation of the curves: dash-dotted CC(stat3), dashed 
CC(ret4), full CC(ret2). Experimental data from M (□), [pi (o), and Isi (•). 
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FIG. 25. Neutron Polarization Py{n) of deuteron photodisintegration 
CC(ret4), full CC(ret2). Experimental data from (•). 
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